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Abstract: We study a two Higgs doublet model extended by vectorlike leptons mixing 
with one family of standard model leptons. Generated flavor violating couplings between 
heavy and light leptons can dramatically alter the decay patterns of heavier Higgs bosons. 
We focus on pp ^ H ^ W where z /4 is a new neutral lepton, and study 

possible effects of this process on the measurements of pp WW and H WW since it 
leads to the same final states. We discuss predictions for contributions to pp WW and 
H WW and their correlations from the region of the parameter space that satisfies all 
available constraints including precision electroweak observables and from pair production 
of vectorlike leptons. Large contributions, close to current limits, favor small tan (3 region 
of the parameter space. We find that, as a result of adopted cuts in experimental analyses, 
the contribution to pp WW can be an order of magnitude larger than the contribution 
to H ^ WW. Thus, future precise measurements of pp WW will further constrain the 
parameters of the model. In addition, we also consider possible contributions to pp WW 
from the heavy Higgs decays into a new charged lepton 64 {H WpVfj), exotic SM 

Higgs decays, and pair production of vectorlike leptons. 
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1 Introduction 

Among simple extensions of the standard model (SM) are those with extended Higgs sector 
and extra vectorlike leptons near the electroweak (EW) scale. Since masses of vectorlike 
leptons are not related to their Yukawa couplings, in the absence of mixing with SM leptons, 
they are not strongly constrained by experiments. However, even small Yukawa couplings 
between SM leptons and vectorlike leptons can significantly affect a variety of processes 
and can dramatically alter the decay patterns of heavier Higgs bosons. 

We consider an extension of the two Higgs doublet model type-H by vectorlike pairs 
of new leptons: SU(2) doublets SU(2) singlets Er^r and SM singlets Nr^r^ where 

Lr and Er have the same hypercharges as SM leptons. We further assume that the new 
leptons mix only with one family of SM leptons and we consider the mixing with the 
second family as an example. The mixing of new vectorlike leptons with leptons in the 
SM generate flavor violating couplings of VE, Z and Higgs bosons between heavy and light 
leptons. These couplings can result in new decay modes for heavy CP even (or CP odd) 
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Higgs boson: H and H e 4 /i, where 64 and z /4 are the lightest new charged and 

neutral leptons. These decay modes, when kinematically open, can be very important, 
especially when the mass of the heavy Higgs boson is below the tt threshold, about 350 
GeV, and the light Higgs boson (/i) is SM-like so that H ZZ, WW are suppressed or 
not present. In this case, flavor violating decays H ot H ^ compete only with 

H ^ bb (for sufficiently heavy H also with H hh) and can be large or even dominant. 
Subsequent decay modes: 64 ^ 64 ^ Z/i, 64 ^ hji and z /4 ^ Vk/i, z /4 ^ Zz/^, 

z /4 ^ lead to many possible final states. 

In this paper, we focus on pp ^ H and study possible effects of 

this process on the measurements of pp WW and H WW since it leads to the 
same final states. This process was previously studied by us in a model independent 
way [ 1 ] in the connection with the ATLAS excess in pp WW [ 2 ]. The results were 
presented in terms of the Higgs mass, the mass of z /4 and the product of branching ratios 
BR(iL ^ i^A^ij.) BR(z /4 ^ Wp). Here we study the process in detail in the two Higgs 
doublet model type-H which is perhaps the simplest realization of the scenario. We discuss 
predictions for contributions to pp WW and H WW and their correlations from 
the region of the parameter space that satisfies all available constraints including precision 
electroweak observables [3] and constraints from pair production of vectorlike leptons [4]. 
Large contributions, close to current limits, favor small tan /3 region of the parameter 
space. We find that, as a result of adopted cuts in experimental analyses, the contribution 
to pp WW can be an order of magnitude larger than the contribution to iL ^ WW. In 
addition, we also consider possible contributions to pp WW from H cap VL/iz/^, 
from similar processes involving SM-like Higgs boson and from pair production of vectorlike 
leptons. 

Vectorlike leptons near the electroweak scale provide a very rich phenomenology and 
were studied in a variety of contexts. Most of the previous studies would apply also to 
the two Higgs doublet model we consider here since we assume type-H couplings of Higgs 
doublets to fermions relevant for supersymmetric extensions and we also consider the limit 
when the light Higgs is SM-like which is relevant for SM extensions by vectorlike fermions. 
For example, analogous processes involving SM-like Higgs boson decaying into 2£2u or M 
through a new lepton were previously studied in ref. [5] and the 41 case also in ref. [ 6 ]. 
Possible explanation of the muon g-2 anomaly with vectorlike leptons was studied in [7, 8 ]. 
Further extensions with vectorlike quarks and possibly Z' are straightforward and offer 
possibilities to explain anomalies in Z-pole observables [9-12]. In addition, extensions 
with complete vectorlike families were considered that provide an understanding of val¬ 
ues of gauge couplings from IR fixed point behavior and threshold effects of vectorlike 
fermions, as in insensitive unification [13, 14]. Many studies were also done in supersym¬ 
metric framework, see for example Refs. [15-20], and in various frameworks the constraints 
from precision electroweak data have been analyzed [21-27]. Further discussion and more 
references can be found in a recent review [28]. 

This paper is organized as follows. In section 2 we present two Higgs doublet model 
type-H with vectorlike leptons mixing with one family of the SM leptons and derive formulas 
for couplings of Z, W and Higgs bosons to leptons. In section 3 we discuss branching 
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ratios of the heavy Higgs boson H and neutral lepton z /4 and find approximate expressions 
assuming small mixing between heavy and light leptons. The scan over the parameter space 
of the model is described in section 4 together with constraints imposed from precision 
electroweak data and direct searches for new leptons. The approximate formulas derived 
in previous sections are useful for understanding the main results which are presented in 
section 5. For completeness, we also discuss possible contributions to pp WW from 
H e 4 /i ^ Wfi^n^ in section 6 , from the SM-like Higgs boson in section 7 and from 
pair production of vectorlike leptons in section 8 . We summarize and present concluding 
remarks in section 9. 

2 Model 

We consider an extension of a two Higgs doublet model by vectorlike pairs of new leptons: 
SU(2) doublets SU(2) singlets Er^r and SM singlets Nr^r. The quantum numbers 

of new particles are summarized in Table 1. The Lr and Er have the same quantum 
numbers as the muon doublet pR (we use the same label for the charged component as for 
the whole doublet) and the right-handed muon pR respectively. We further assume that the 
new leptons mix only with one family of SM leptons and we consider the mixing with the 
second family as an example. This can be achieved by requiring that the individual lepton 
number is an approximate symmetry (violated only by light neutrino masses). The results 
for mixing with the first or the third family could be obtained in the same way. The mixing 
of new leptons with more than one SM family simultaneously is strongly constrained by 
various lepton flavor violating processes and we will not pursue this direction here. Finally, 
we assume that leptons couple to the two Higgs doublets as in the type-H model, namely 
the down sector couples to Hd and the up sector couples to Hu- This can be achieved by 
the Z 2 symmetry specified in Table 1. The generalization to the whole vectorlike family of 
new leptons, including the quark sector, would be straightforward. 

Table 1. Quantum numbers of standard model leptons, extra vectorlike leptons and the two Higgs 
doublets. The electric charge is given by Q = T 3 + T, where T 3 is the weak isospin, which is +1/2 
for the first component of a doublet and - 1/2 for the second component. 



Ml 


Ll,r 

El,r 

Nl,r 

Hd 

Hu 

SU(2)l 

2 

1 

2 

1 

1 

2 

2 

U(1)y 

1 

2 

-1 

1 

2 

-1 

0 

1 

2 

1 

2 

Z 2 

+ 

- 

+ 

- 

+ 

- 

+ 


With these assumptions, the most general renormalizable Lagrangian containing Yukawa 
and mass terms for the second generation of SM leptons and new vectorlike leptons is given 
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£ D - y^flL^RHd - XEflLERHd - XLLLynHd - XLeErH^ - XH^^ErLr 

- rnP^lErHu - rLrNrHu - RH^NrLr (2.1) 

— MlLrLr — MrErEr — MnNrNr + h.c. , 


where the first term is the usual SM Yukawa coupling of the muon, followed by Yukawa 
couplings to Hd (denoted by various As) that will result in masses and couplings of the 
charged leptons, Yukawa couplings to (denoted by various rs) that will result in masses 
and couplings of the neutral leptons, and finally mass terms for vectorlike leptons. The 
components of doublets are labeled as follows: 





Ll,r 





( 2 . 2 ) 


where the neutral Higgs components develop the vacuum expectation values (H^) = Vu 
and (H^) — Vd- We assume that both are real and positive as in the CP conserving two 

Higgs doublet model with ^JvuE Vd = v = 174 GeV and we define tan/3 = Vulvd- 

After spontaneous symmetry breaking the resulting mass matrices in the charged and 
neutral sectors can be diagonalized and we label the two new charged and neutral mass 
eigenstates by 64 and 65 and z /4 and z /5 respectively. Couplings off all involved particles 
to the Z, W and Higgs bosons are in general modified because SU(2) singlets mix with 
SU(2) doublets. The flavor conserving couplings receive corrections and flavor violating 
couplings between the muon (or muon neutrino) and heavy leptons are generated. The 
couplings resulting from the mixing in the charged sector were discussed in detail in ref. [ 8 ] 
in the connection with the muon g-2 anomaly. Here we will focus on couplings resulting 
from the mixing in the neutral sector. These are also more relevant for the discussion of 
the contribution of the Higgs boson decays to pp WW. 

The mass matrix in the neutral lepton sector is given by: 


(i>M Ll Nl)m, 


^ Ur = 0^ 

L% 

Nr J 


(^.L^lNl) 


^ 0 0 KjvVu ^ 

0 Ml kVu 

yo KVu Mn y 


l^R = 0\ 

L% 

Nr J 


(2.3) 


where we inserted er = 0 for the right-handed neutrino which is absent in our framework 
in order to keep the mass matrix 3 x 3 in complete analogy with the charged sector. For 
the discussion of couplings it is convenient to define vectors = (^/x, A^y,)^ and 

i/Ra = (er = 0, Nr)^. The mass matrix Mjy can be diagonalized by a biunitary 

transformation 


vIm.Vr 


/o 0 0 \ 

0 0 
y 0 0 y 


(2.4) 


resulting in masses for z /4 and z /5 leaving the muon neutrino massless. The light neutrino 
masses can be generated by a variety of ways. Once they are generated, the mixing of 
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light neutrinos with vectorlike leptons results in corrections to both the masses and mixing 
angles controlled by Yukawa couplings in eq. (2.1). 

For better understanding of corrections to gauge and Yukawa couplings discussed later, 
approximate analytic formulas for diagonalization matrices are useful. These can be ob¬ 
tained in analogy with those in the charged lepton sector given in ref. [8]. In the limit 

KJV^U, KVu, KVu < Ml, Mn (2.5) 


with Ml and Mjq not close to each other, we find 


and 


Vl 


1 _ 

2M% 

MlMn 

_ I^N^u 

Mn 


K,NVu kMl+kMn 
Ml M%-Ml 
-I _ {MlR^Mni^)‘^Vu 
2{M%-MlY 
{Mli^^Mnk)vu 
Ml-Ml 


KnVu \ 

Mm 

{MlK+Mmk)Vu 
Ml-Ml 

, _ KN^vj _ {Mlk+MmkYvI 
2Ml 2{Ml-MlY j 


Vr 


/l 0 

n 1 _ {Mlk+MmkYvI, 
2{Ml-MlY 
^ {MlI^^MnI^)Vu 

Ml-Ml 


0 \ 

{Mlk,+Mnk,)vu 

Ml-Ml 

1 _ {Mlk+Mmk,YvI 
2{Ml-MlY / 


( 2 . 6 ) 


(2.7) 


up to corrections of 0{e^) where e = {kn, K,ii)vu/{ML, M^). The mass eigenvalues are 
0, ML + 0{e^), Mn + 0{e^). However, in our numerical analysis we do not use any approx- 
imations. 


2.1 Couplings of the Z and W bosons 

Couplings of the muon and new heavy leptons to the Z and W bosons are modihed from 
their SM values because SU(2) singlets mix with SU(2) doublets. These couplings can be 
written in terms of Vl and Vr^ dehned in eq. (2.4), and of the analogue matrices Ur and 
Ur that are related to the charged lepton sector and that were discussed in detail in ref. [8] 
(with the replacement v ^ due to the two Higgs doublet model). The couplings of the 
Z boson to charged leptons can be found in ref. [8] and those to neutral leptons follow from 
the kinetic terms: 


^kin ^ T^La0a^La + T^Ra0a^Ra 

= ^La{Vl)acilPc{^L)chi^Lh + (H^)ac0c(^i^)c6^i^65 (2.8) 


where the vectors of mass eigenstates are URa = (^/x, ^L 45 ^Ls)^ and similarly for URa = 
(Or = 0, i>i^ 5 )^. We label the components of vectors and diagonalization matrices by 

2, 4 and 5 because they correspond to 2nd, 4th and 5th mass eigenstate. The covariant 
derivative is given by: 


D 


fia 


— f) — j _ i. _ T^Z 

^ cos 0 m “ 


(2.9) 
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where the weak isospin is +1/2 for neutral components of SU(2) doublets and 0 for 
singlets. Defining couplings of the Z boson to leptons fa and as 

c D , ( 2 . 10 ) 

we find: 


ZUa^b _ 

Ol — 

Zya^h _ 

Or 


9 


2 cos 9\y - 


iVl)a2iVL)2b + {Vl)a4{VLU 


g 


(V^)a4(VR)46 . 


2 cos 0\Y 

The usual SM couplings of left-handed neutrinos, 

Zya^h^ 9 


{gr^^)sM = 


2 cos 6w 


ab 


are modified by 


Sg 


ZVaVb 


g 


{Vl)ab{VLU. 


“ 2 COS 

The couplings of the W boson originate from the kinetic terms: 

Ckin D ^ + Lh^Ll + L-^) W;^ + h.c. 

= + (l^La{Vl)a27'^{UL)2b^Lb + l^La{Vl)a47'^{UL)4beLb 


+ 


l^Ra{V^)a4'y'^{UR)4beRb'^ 1+ + h.C. , 


( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


where clb, (eRb) are the charged left-handed (right-handed) mass eigenstate and Ur, Ur 
are the corresponding diagonalization matrices [5]. Defining conplings of the W boson to 
neutrinos Ua and charged leptons ei as 

^ {i>Lal^gY’'’'‘'’’^Lb + i>Ra'y'^gR’'‘"''^eRb^ W+ -b h.c., (2.16) 

we find: 


WyaCb 

9l 


WyaCb 

9r 


)a2{UL)2b+{Vl)a4{UL)4b 


(2.17) 

(2.18) 


2.2 Couplings of the Higgs bosons 

As a consequence of explicit mass terms for vectorlike leptons, the usual relations between 
the mass of a particle and its coupling to Higgs bosons do not apply. The couplings of 
neutral Higgs bosons to neutral leptons can be obtained from the following Yukawa terms 
in the Lagrangian (2.1): 


D b'j^a Yyab ^Rb h.C. 

— ~^La{V^)acYjycd {^R)db ^Rb + h.C.^ 


(2.19) 
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where 




^0 0 
0 0 K. 

^0 ^ 0 y 


( 2 . 20 ) 


We assume a CP conserving two Higgs doublet model in the limit with the light Higgs h 
being fully standard model like in its couplings to gauge bosons and the heavy CP even 
Higgs H having no couplings to gauge bosons. The mass eigenstates h and H in this limit 
are related to doublet components as follows (see for example ref. [29]): 


/ /i \ / cos/3 sin/3\ / y/2(ReH^ — Vd)\ 

Y —H j — sin /3 cos y y ~ ^w) j 


The Yjy matrix is not proportional to the mass matrix given in eq. (2.3) and thus the Higgs 
couplings are in general flavor violating. Defining couplings of mass eigenstate leptons fa 
and fi) to CP-even Higgs bosons by 

-C D /La fmh-^ fLa f Rb H + h.c., (2.22) 

we find: 


Xt^,^=smf3iVlY,VR)ab, 

-Af^,^=COS/3(V/VV?)a6. 


Since Yj^Vu = Mjy — diag{0^ Ml^ the Higgs boson couplings can be also written as: 



/ 0 

0 

0 

\ 

II 

0 


0 



(0 

0 

m^5 

/ 


/ 0 

0 

0 

\ 


0 

rrii,^ 

0 



(0 

0 


/ 




/o 0 0 \ 

0 Ml 0 
^0 0 Tf/v j 

/o 0 0 \ 

0 Ml 0 

^0 0 Tf/Y y 




Vr, 


(2.24) 


where we used Vu — vsin/S. The first terms in above equations represent the expected 
relations between fermion masses and their couplings to Higgs bosons and the second term 
represents contributions from the terms. This form of couplings makes it obvious 

that in the absence of vectorlike masses the couplings of h to leptons are fully SM-like, 
while couplings of H are enhanced by tan /3 as expected in the limit we assume. 

Couplings to charged leptons follow from terms in eq. ( 2 . 1 ) and can be obtained 
from those in eqs. (2.23) with replacements: Vl^r W ^ Ye /3 ^ /3 + 7 r/ 2 , see 

also ref. [ 8 ] in the case of SM. The corresponding formulas to eqs. (2.24) would show that 
couplings of h have the usual SM strength, up to contribution from while couplings 

of H to charged leptons are suppressed by tan^S. Finally couplings of the CP-odd Higgs 
boson, A, copy those of H up to the usual 75 factor. 
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3 Branching Ratios 


We collect expressions for the relevant branching ratios for the process pp ^ H ^ 

W and provide several approximate formulas in the limit of small mixing between neu¬ 
tral leptons discussed in the previous section. These formulas will be useful for qualitative 
understanding of results. From now on, we drop the hat notation for mass eigenstates and 
also label the mass eigenstates i)l 2 and 62 as and p. 

Sizable decay modes of the heavy CP even Higgs boson are 2 ^- 4 bb and gg for mn < 
250 GeV. As discussed in the previous section we assume that H does not have direct 
couplings to pairs of gauge bosons and that decay modes to other Higgs bosons are not 
kinematically possible. However, our results could be straightforwardly modified to account 
for additional sizable decay modes of H. 

The partial decay width of H ^ (where we include both z/ 4 Z/^ and z/ 4 Z/^ final 

states) is given by: 

(3,1) 

where 


^^ 1^4 = -^(VjJ)24(Vr)44 + —^{Vl)2h{VR)u 

z^w'^sin/3cos/3 |~_ ^ -\-MjsiR) -\-M]s[RYv‘^ siv? (3 

~ jPi Ml-Mf 2(Ml - M2)2 


(3.2) 
, (3.3) 


where the second line is an approximate formula in the limit of small mixing discussed in 
section 2. Note, that this limit assume the z /4 is mostly the doublet with mass originating 
from Ml. For an approximate formula corresponding to a singlet-like neutral lepton, the 
should be used instead. This coupling is given by 


A 


H 




cot 13 


V 


— KjSf COS /3 


(H2)24(Pr)45 H-^(H2)25(Hr)55 

{Mli^ -h MnR)‘^v‘^ sin^ /3 Rv‘^ sin^ /3 Mli^ + M^R 
2 (M2^ - MlY Mtv M%-Ml \ ’ 


(3.4) 

(3.5) 


where the second line is an appropriate approximate formula in the case of singlet-like 
lepton with mass originating from Mw- 

The decay width of H ^ bb is given by 


V{H bb) 


3Gf 

4:^/27T 


mHml{mH) tan^ ^ [l + \q + A|^] , 


(3.6) 


where mAmn) is the running b-quark mass evaluated at the scale mn and the correction 
factors Agg and A|^ can be found in ref. [30] . The decay width of H ^ gg is given by 


r{H ^ gg) 


Gpa^mL cot^ (3 
36\/27r^ 
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3 

4 


^ 1/2 


2 


(3.7) 






















with 


^ 1/2 — 2 


r + (r - l)/(r) 


(3.8) 


where r = m^/4m^ and /(r) = arcsin^ foi" t <1. 
The branching ratio of ^ z/ 4 Z/^ is then given by 


BR{H = 


V{H 1 ^ 4 ^ 12 ) 


r{H ^ U4iy^) + r{H ^ 66 ) + r{H ^ 


(3.9) 


The neutral lepton z /4 can decay into standard model leptons and the Higgs, IT, and 
Z bosons. Neglecting the muon mass, the partial decay width of z /4 ^ is given by: 


r(z/4 ^ hi^n) - (^1 


m 


2 \ 2 


m 




where 


^- ■^{VI)2a{Vr)aA - ■^{Vl)2b{VR)^A- 


The partial decay width of z /4 ^ IT/i is given by: 


r(z24 ^ hF+/i“) = 


m 


UA 


327r 


{ 9 . 


WuAlJi \2 


+ { 9 ^ 


Wl'AlJi \2 ^VA 


J- o 


m 


UA 


1 + 2 - 


Ml 


w 


m 


UA 


where 


{VI)a2{Ul)22 + {VI)aa{Ul)a2 


WvAtx _ _9_ 

V 2 L 

WUA, _ . 




9r - ^y'R> 

Finally, the partial decay width of z /4 ^ is given by: 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


r(i/4 ^ Zv.) = L 

\ fA) 32^vyL ; 


M\ 


2\2 


m 


UA 


1 + 2 - 




m 


UA 


where 


ZlZa^A 

9l 


9 


{VI)a2{Vl)22 + {VI)aa{Vl)a2 




(3.15) 


(3.16) 


2 cos 9y[z 

Assuming only these decay modes of z/ 4 , the branching ratio of z /4 ^ VF/i is given by 

r(z/4 ^ iF+/i“) 


BR(l/4^hh+M ) Y{RA^hR^)+T{Ri^W+^JL-) + T{RA^ Zu^) 
The branching ratio of ^ VFis defined as 

BR{H WpLUf,) = BR{H BR(i/4 ^ TR+zi"). 


(3.17) 


(3.18) 
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4 Parameter space scan and constraints from precision electroweak data 
and direct production 


We perform a scan over all the model parameters introduced in section 2 over the ranges 


Ml^n e [0,500]GeV , (4.1) 

G [—0.5,0.5] , (4.2) 

tan/3 G [0.3, 3] . (4.3) 


We fix the mass term of the SU{2) singlet charged vectorlike lepton Me — 1000 GeV. 
We simplify the decay patterns of the heavy Higgs by requiring > uih (to avoid 
H Z/ 5 X channels) and rrijy^ > mHl‘2 (to avoid decays into pairs of heavy vectorlike 
leptons). Moreover we include mixing exclusively in the neutral sector. 

We impose constraints from precision EW data related to the muon and muon neutrino 
that include the Z pole observables (Z partial width to the invisible width, forward- 

backward asymmetry, left-right asymmetry), the W partial width, and the muon lifetime. 
We also impose constraints from oblique corrections, namely from S and T parameters. 
Unless specified otherwise these are obtained from ref. [3]. Finally, we impose limits from 
direct searches: the LEP limits on masses of new charged leptons, 105 GeV, and the limits 
on pair production of vectorlike leptons at the LHC summarized in ref. [4]. Constraints 
on the production of heavy Higgs will be discussed in the following section together with 
results. 

Constraints on the muon couplings were already discussed in ref. [ 8 ]. Precision elec¬ 
troweak measurements constrain modification of couplings of the muon to the Z and W 
bosons at ^ 0 . 1 % level which, in the limit of small mixing, approximately translates into 
95% C.L. bounds on Xe^l couplings: 


XEVd 

Me 


< 0.03, 


XLVd 

Ml 


< 0.04, 


(4.4) 


assuming only mixing (Yukawa couplings) in the charged sector. 

In the neutral lepton sector the strongest limits are obtained from the muon lifetime. 
In what follows we discuss this limit together with the invisible widths of the Z boson and 
constraints from direct production of vectorlike leptons. 


4.1 The muon lifetime 


The Fermi constant Ge determined with a high precision from the measurement of muon 
lifetime. In the standard model Ge = {y/2l8)g‘^ jM^ while in our model one of the g/V2 
factors is replaced by gY^^^ given by 




9 

V2 


{VI)22{Ul) 22 + iVl)24iUL)42 


(4.5) 


The allowed range for is obtained from the uncertainty in the W mass, Mw = 

80.385 ± 0.015 GeV. The relative uncertainty in is 2 x = 3.73 x 10“^ and we set 
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the 95% C.L. upper limit on the deviation of from gl\f2 as: 


9l 


9l^ 


- 1 


< 2 X 3.73 X 10“^ = 7.46 x 10“^ . 


(4.6) 


Assuming no mixing in the charged sector and using the small mixing approximation 
eq. ( 2 . 6 ), 


wuu,ij. _ g 




„2 

1% Ar U, 


and we obtain an approximate 95% C.L. upper bound on the size of kjv conpling: 

KnVu 


Mat 


< 0.04. 


(4.7) 


(4.8) 


Considering also mixing in the charged sector, the bound is shared between hcn and A^: 


The partial decay width of VC ^ /j.u^ depends quadratically on the coupling. 

However, it is measured with about 2% precision and thus the resulting constraint on the 
coupling is significantly weaker. 


4.2 Invisible width of Z 


The partial width of Z ^ is given by 


where 


r(z ^ 


247r^^^ 


Zu^u^ _ 

9l — 


9 

2 cos 9w 
9 

2 cos 9w 


(V2)22(Hl)22 + (H2)24(Hl)42 


1 - 


Ml 


(4.10) 


(4.11) 

(4.12) 


where the second line is an appropriate approximate formula in the case of small mixing. 
In this limit, the upper bound on obtained from muon lifetime, eq. (4.8), suggests 
that r(Z ^ can be modified at most at 0.3% level (the sign of the correction 

is always negative). Since we assume that only one generation of SM leptons mix with 
vectorlike pairs, the invisible width of Z can be lowered at most by 0.1%. This is also 
visible in figure 1 where we consider randomly generated points in the ^5 ^5 Ml and 
Mat parameter space for fixed mn = 155 GeV and rriiy^ = 135 GeV (different choices 
of masses do not sizably affect the allowed ranges), assuming no mixing in the charged 
lepton sector, and impose all EW precision constraints (including direct productions bounds 
discussed in section 4.3 below). In the left and right panels of figure 1 we consider the 
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— 0.020 


0.015 



0.000 0.005 0.010 0.015 0.020 


Figure 1. Parameter space scan points (for mn = 155 GeV and rriy^ = 135 GeV) that survive all 
EW precision constraints. The blue, cyan, magenta and red points have singlet fraction (defined as 
|(^L)45p/^ + I(^A)45|^/^) ranges [0,5]%, [5,50]%, [50,95]%, and [95,100]%, respectively. In 

the left panel, we show the ~ I plane; on the right axis we show the 

corresponding ratio. In the right panel we show the off-diagonal and 

gauge couplings. 


and -bf ^"1 planes. 

Both the upper limit on tv]\fVu/M]\f^ eq. (4.8), and the resulting largest possible effect in 
r^inv follow closely those obtained from the approximate formulas. Points with z /4 being 
mostly singlet (red points) or mostly doublet (blue points) cluster very near the line that 
assumes the approximate relation in eq. (4.7) is exact and even highly mixed scenarios 
(cyan and magenta) are not very far. 

Since the ratio of the measured value of the Z-boson invisible decay width and its SM 
expectation is 



499.0 ±1.5 (MeV) 
501.66 ±0.05 (MeV) 


0.995 ± 0.003 , 


(4.13) 


the invisible width does not provide additional constraint to that obtained from the muon 
lifetime. This however assumes that only the 2nd generation of SM leptons mixes with 
vectorlike leptons. The effect on invisible width can be larger if more generations mix with 
vectorlike leptons or if one considers mixing with the 3rd generation instead of the 2nd 
since the constraints on 3rd generation couplings are weaker. 

In conclusion, couplings of SM gauge bosons to the second family of leptons, 
and ^ can deviate from their SM values by less than ^ 0.1%. Moreover, within the 

explicit model we consider, these constraints imply upper limits of order ^ 0.02 on the new 
off-diagonal gY^^^ gauge couplings as we can see in the right panel of figure 1 . 


4.3 Direct production of vectorlike leptons 

Let us first consider constraints from Drell-Yan production of Z/ 4 Z /4 or z/ 4 e 4 leading to at 
least 3 leptons in the final state and The cross section of pp Z/ 4 Z /4 is proportional 
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[GeV] 

105 

125 

150 

200 

300 

^ BR (z/4 y lR)\max 

0.090 

0.141 

0.141 

0.164 

0.582 

\Re4V4^ BR(z /4 y g^\max 

0.109 

0.203 

0.267 

0.355 

1 


Table 2 . Upper bounds on BR^(z/4 Wji) and Re^y^x BR(z/4 Wji) obtained from 

ref. [ 4 ] for several masses of z/4. The limits on R^^y^x BR(z/4 Wfi) assume that BR(e4 ^ 
= 1 and rrie^ = rriy^. 



Figure 2 . Upper bounds on Ry^^x BR(z/4 Whu^) as functions of rriy^. 


to where the couplings are given in eqs. (2.11)-(2.12). Thus the cross 

section is modified from the one that corresponds to fully doublet z/4 by factor 


Ru^ua. — 2 


(^i)42(Vz.)24 + (^i)44(Vz.)44 


2 1 

+ - 


(V^)44(Vr)44 


1 2 


(4.14) 


At present, searches for anomalous production of multilepton events constrain only the 
case when both Z/4S decay to VU/i and the limits on Ry^y^x BR^(z/4 Wfi) can be read 
from Table 2 of ref. [4]. They are summarized in our Table 2 for reference values of rriy^. 
In our numerical analysis we interpolate these results for other values of rriy^. 

Similarly, the cross section of pp e4Z/4 is proportional to where 

the couplings are given in eqs. (2.17)-(2.18). Thus the cross section is modified from the 
one that corresponds to fully doublet z/4 and 64 by factor 

^ Vi)42(eL)24 + (U)44 


R. 


'€41^4 — 


(Ul) 


44 


2 1 
+ 2 


{V^)u{Ur)m 


(4.15) 


The analysis in ref. [4] , assuming rrie^^ = rriy ^, shows strong limits on this production cross 
section. Some decay modes of 64 and z/4 are consistent with data only for masses higher 
than 500 GeV. In our analysis we are focusing on the case with no mixing in the charged 
lepton sector. In this case BR(e4 ^ Wix^) — I (the required coupling originates from 
mixing in the neutral sector) and the limits on R^^y^x BR(z/4 Wfi) can be obtained 
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from Table 2 of [4]. Other decay modes of z/4 are not constrained in this case. The upper 
bounds are summarized in Table 2. We again interpolate these results for other values of 


Finally, let us comment on a single production of a new lepton. The W — U 4 — jj. 
coupling results in a production of z/4 through the process pp VF* ^ z/4/i which also 
can lead to at least 3 leptons with in the final state. The bounds were discussed in 

ref. [31] in the context of TeV scale seesaw models with very small lepton number violating 
terms, see for example ref. [32], which are not constrained from the same-sign dilepton 
searches. 


The cross section of pp z/4/i is proportional to where the 


cou¬ 


plings are given in eqs. (2.17)-(2.18). Thus the cross section is modified from the one that 
corresponds to full strength coupling of the two leptons to W by factor 




(vl 




(Ul) 


22 


+(v^) 


lM4 


(Ul) 


12 


42 


+ 


{Vi)44{UR)42 


12 


(4.16) 


We closely follow ref. [4] to set limits from the ATLAS searches for anomalous production 
of multilepton events [33] on three decay modes: 


pp VF* ^ z/4/i ^ VF/i/i , , hu/^p , (4-17) 

where h is the SM Higgs with mass 125 GeV. The obtained upper bounds on 
BR(z /4 ^ Wp^ are shown in Fig. 2 as functions of We see that for any 

combination of branching ratios the constraint on is at most of 0(10“^). This limit is 
much weaker than those obtained from precision EW data; in fact, for the surviving points 
in figure 1 the maximum value of Ry^^ is of 0(10“^). 


5 Main results: contributions oi H pv to pp ^ FFFF and H FFFF 


In this section we explore the impact that this model has on pp (VFVF, H WW) 
iul'Ty' measurements. We show detailed results for the two representative points (m/f, rriy^) = 
(155 GeV, 135 GeV) and (250 GeV, 230 GeV) in figures. 3 and 4. In figures. 5-7 we show 
how these results vary for different values of mn and rriy ^. 

In figure 3 we present the results of the scan described in section 4 for the reference point 
= (155 GeV, 135 GeV) discussed in ref. [1]. The blue, cyan, magenta and red 
points have z/4 with singlet fraction (1(^2)451^/^+l(^^)45l^/2) ranges [0,5]%, [5,50]%, 

[50,95]%, and [95,100]%, respectively (note that in some of these plots blue/cyan/magenta 
colors are not easily distinguishable). In the two upper plots the black contours are the 
values of the effective pp WW cross section as defined in [1] for the epu^i'i^ ([^Np^]e/i) 
and pph'^h '12 final states, respectively. In parenthesis we show the corresponding 

effective pp ^ H ^ WW cross sections • These effective cross sections ^ 

are explicitly defined as: 


(JNP 


(j{pp H ^ W£jy ijyi'jy') Anp 
p BR{W ivY Asm 


^An extended discussion of the effective cross sections is presented in section 2 of ref. [1]. 


(5.1) 
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where r] = 2(1) for e/i (/i/i) final states and the NP and SM acceptances and Asm 
are calculated using the experimental WW and H WW cuts (for the latter we follow 
ref. [34] and consider the six Higgs mass hypotheses discussed in ref. [35] and show the most 
constraining effective cross section). Note that points displayed in the two upper panels 
are identical and that the only difference lies in the contours that depend crucially 

on the very different acceptances for e/i and /i/i final states as well as the factor rj. Note 
that eq. (5.1) implies 


H^WW 

^NP 


aU aWW 
^NP ^SM 

An aww 
^SM ^NP 


WW 

^NP 


(5.2) 


The product of acceptances in this equation is the crucial parameter that controls the size 
of contributions to pp WW that are allowed hy H ^ WW searches. For most (but 
not all) masses that we consider, this ratio is of order 10% (typically A^^p/Ay^ 0(0.1) 
and 0(1)). The smallness of reason for which we can find 

large cross sections while simultaneously surviving H WW bounds. When the 

difference ttih — is large and rriy^ is small, and ttit increase while mii decreases 

implying a larger ^npM^ ratio. For instance, for mn — 250 GeV and rriu^ = 135 GeV 
we find that this ratio can be as large as 2. 

The yellow shaded area is excluded hy H ^ WW searches. The upper bound on the 
effective cross section is independent of and is given by 


H^WW 

^NP 


< min 

n 


1 


1 

rj BR(1F ^ Ipf ’ 


(5.3) 


where is defined in appendix A of ref. [1]. The measurement of the pp WW cross 
section is very sensitive to NNLO QCD corrections which have not been fully implemented 
in the experimental analysis yet. Following, for instance, the discussion around eq. (1.7) of 
ref. [1], the deviation of the pp WW cross section with respect to the SM expectation 
found by ATLAS [2] and CMS [36] are: 


= (l2.7^Ji) pb = (-0.1 ± 5.3) pb 

I pb I = ('‘■5 ± 7.0) pb 

Since these two results adopt different theoretical setups, we refrain from combining them 
into a weighted average. For this reason do not use pp WW data to constrain our model 
and simply quote the allowed values. 

A prominent feature of figure 3 is that for a doublet-like P4 the product of branching 
ratio BR(i7 ^ W— BR(i7 ^ P4P^) x BR(z /4 Wp) is constrained to be very 

small. This is mainly due to bounds from the multilepton plus searches in the Drell- 

Yan pair production process pp Z/4P4 ^ W^W~p'^p~. This can be understood by 
looking at the middle-left panel of figure 3 where we consider the — BR(p 4 Wp) 

plane. The quantity is defined in eq. (4.14). Here the light colored points are 

obtained without imposing any of the constraints discussed in section 4 and the darker 
colored points are those that survive after imposing the muon lifetime bound. Additional 
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Figure 3. Parameter space scan for ttih = 155 GeV and = 135 GeV. The blue, cyan, 
magenta and red points have singlet fraction in the ranges [0,5]%, [5,50]%, [50,95]%, and [95,100]%, 
respectively. In the two upper plots all constraints are imposed and we focus on the BR(i^ ^ 
Wiiv^) — tan/3 plane for the and final states, respectively. The black contours 

are the values of the effective pp WW and pp ^ H ^ WW cross sections (in pb) defined in 
eq. (5.1). The yellow shaded area is excluded by ^ WW searches. In the middle-left plot, we 
consider the — BR(z /4 Wp). Here the light-shaded points do not satisfy the muon lifetime 

constraint and the impact of multilepton + searches from Drell-Yan pair production process 

pp z^ 4 Z /4 ^ W^W~p~ is indicated by the black curve. In the middle-right plot we show the 
BR(z /4 ^ — BR(z /4 Wp). Here the gray points are excluded by multilepton searches. In 

the two lower plots we consider the BR(H’ ^ z/4Z/^) — BR(z/4 Wp) and planes. 


- 16 - 

























0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 





0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 





Figure 4. Parameter space scan for ttih = 250 GeV and = 230 GeV. See the caption in 
figure 3 for further details. 


constraints from oblique corrections are very strong (especially from the S parameter) but 
in the [BR(i7 ^ VP/iz/^),tan/3] plane they do not modify significantly the overall allowed 
region. 

Bounds from multilepton searches exclude the region above the black contour separat¬ 
ing the surviving points in two disconnected regions at low ^^74^/4 with BR(z/4 Wfi) ^ 70% 
and large ^ 1 with low BR(z/4 ^ W 11 )? 

^Note that these arguments rely strongly on the particular choice of masses {{mH— 
(155 GeV, 135 GeV) in this case); a completely different situation characterizes the configuration presented 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 


BRiH^Wnu ^) 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 


BRiH^Wfiu ^) 


Figure 5. Projection onto the [BR(i7 ^ tan/3] plane of the parameter scan described in 

the main text for various values of ttih and . The black contours are the values of the effective 
pp WW and pp ^ H ^ WW cross sections (in pb) for the final state. See the caption 

in figure 3 for further details. 


At small i?iy4i/4 the z/4 is mostly singlet, the second term in eq. (3.16) is suppressed by a 
factor (V^)44 with respect to the first and the 1 ^ 4 , —Z — coupling is controlled by the single 
quantity (V^)42 (we remind the reader that (Vl )22 is very close to 1). Under the assumption 
of no mixing in the charged sector, the matrix U is the identity and the — W — fi coupling 
in eq. (3.13) is also controlled by the parameter (u2)42- As a consequence the ratio of these 
two couplings is the same as in the SM (i.e. independent of flavor mixing parameters), 
implying an almost constant z/4 ^ VP/i branching ratio (^ 70%). At large the z/4 is 

mostly doublet, both terms in the — Z — coupling are of similar size, and the ^ 4 ^ pi 
branching ratio can acquire any value depending on the choice of input parameters. On 
top of this one should note that the z/4 ^ channel is phase space suppressed for the 
case rriy^ — 135 GeV. These considerations are also illustrated in the middle-right plot of 
figure 3 where we show the points in the BR(z/4 ^ Zi^fj^) — BR(z/4 -^Wpi) plane. Here the 
gray points are excluded by multilepton searches and, to a lesser extent, oblique corrections. 

The surviving region at large is also characterized by a very small H — ^ 4 ^ — 

in Fig 4 and discussed later on. 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 


BRiH^Wfiu ^) 


Figure 6. Projection onto the [BR(i7 ^ tan/3] plane of the parameter scan described in 

the main text for various values of ttih and . The black contours are the values of the effective 
pp WW and pp ^ H ^ WW cross sections (in pb) for the final state. See the caption 

in figure 3 for further details. 


coupling as we can see in the lower-left panel of figure 3. In fact, an almost completely 
doublet z/4 requires very small couplings k and implying a strong suppression of the 
Yukawa coupling given, for doublet z/4, in eq. (3.3). This can be seen in the lower-right 
panel of figure 3 where we show the values of the Yukawa couplings A^^^ and for the 
points that survive all constraints. Therefore BR(i7 ^ ^ and hence BR(i7 ^ WpLU^)^ 

are very small for doublet-like z/4. If z/4 is singlet-like, the SM-like Higgs Yukawa coupling 
is given, in the limit of small mixing, by A[5^^^ ^ A^wsin/3, see eq. (3.5). In this case, at 
fixed Mw, the muon lifetime limit (4.8) translates into a direct constraint on the Yukawa 
coupling 

In figure 4 we present the = (250 GeV,230 GeV) case. Now, the large z/4 

mass implies that the decay mode z/4 ^ is no more phase space suppressed and can be 
dominant in large part of the parameter space as we can seen directly in the middle-right 
plot in figure 4 and indirectly in the middle-left plot where BR(z/4 -^Wfi) can only be as 
large as 60%. On top of this, the constraint from multilepton -h searches is weaker 

(this happens generally for rriy^ > 150 GeV as we can see in Table 2), implying that there 
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Figure 7. Projection onto the [BR(z/4 ^ Zz/^),BR(z/4 ^ W jj.)] plane 

of the parameter scan described in the main text for {TnH^TUjy^) = 

(140 GeV,135 GeV),(250 GeV, 135 GeV),(155 GeV, 125 GeV) and (155 GeV, 150 GeV). 
See the caption in figure 3 for further details. 




Figure 8. In the left plot we show how the envelope of the points changes for fixed rriy^ = 135 GeV 
and rriH ^ [140, 250] GeV. In the right plot we take mn = 155 GeV and rrij^^ G [125,150] GeV. 
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is a large region of allowed parameter space in which the z/4 is mostly doublet as can be 
seen in the two top plots in figure 4. Even though there are many points for which the z/4 
doublet fraction is large, the corresponding values for BR(i7 ^ W fiv) are much smaller 
than for typical singlet points. This is because the H — 1^4 — coupling for a doublet-like 
z/4 is suppressed compared to the singlet-like z/4 by RVu/Mn^ see eqs. (3.3) and (3.5). From 
the bottom-right plot in figure 4 we see that the actual bounds on and are 0.05 
and 0.17, respectively. Given that the ratio of these couplings is equal to tan/3, the second 
bound is effectively set by the perturbativity request tan/3 > 0.3. 

In figures 5 (for the e/i final state) and 6 (for the fifi final state) we present the result of 
similar scans for (m/f, rrijy^) — (140 GeV, 135 GeV), (250 GeV, 135 GeV), (155 GeV, 125 GeV) 
and (155 GeV, 150 GeV). The interpretation of these plots is similar to that of figure 3. The 
main difference between these plots is the maximum value allowed for BR(i7 ^ 

In figure 7 we show the BR(z/4 ^ — BR(z/4 -^Wfi) plane for each set of masses. 

In figure 8 we show the envelopes of the allowed parameter space for a wide range of 
masses; in the left plot we take rriy^ — 135 GeV and uih G [140, 250] GeV and in the right 
plot we have ttih — 155 GeV and rrijy^ G [125,150] GeV. This effect is due to change in the 
phase space available for H W as the masses vary. 

Assuming that the acceptances ratios remain constant 

when increasing the center of mass energy from 8 to 13 TeV, the contours in figures. 3- 
6 will simply scale with pp ^ H production cross section. For instance, for mn = 250 GeV 
the rescaling factor is about 2.776 [37]. 

Finally let us comment on the reach of the next LHC run at 13 TeV with a luminosity 
C = 100 fb“^. Taking into account that a{pp WWY^rj.^Y/a{pp l^l^)8TeV — ^ 
and that the uncertainty on is ^ 5 pb (see eq. (5.4)), we estimate ^ 3 pb. 

Moreover, our new physics contributions to scale with the pp ^ H cross section and 

increase by a factor ^ 2.5 [37] at 13 TeV. Taking these considerations into account, direct 
inspection of figures. 3-6 shows that most of the presently allowed parameter space will be 
tested. For instance, with respect to the top-left panel of figure 3, LHC8 with 20 fb“^ is 
sensitive to points below the 5 pb contour while LHC13 with 100 fb“^ will be sensitive to 
points roughly below the 1.2 pb one (that will correspond to [crJJ^^]i3TeV — 3 pb). 

6 Contributions from H 

In this section we discuss contributions to pp ii'iy£iy£f stemming from heavy Higgs pro¬ 
duction and decay into a charged vectorlike lepton and a muon: 

pp ^ H ^ Cap WVfjiP . (6-1) 

We begin our analysis with a model independent study of this channel along the lines of the 
analysis presented in ref. [1]. Our main results are summarized for the ep and pp modes 
in the two panels of figure 9. These figures are very similar to figure 1 of ref. [1]. The 
blue contours are the values of the effective WW cross section that we obtain for 

BR(i7 ^ Wpi'^)cot^ P = 1 (in this section only we define BR(i7 ^ Wpv^) = BR(i7 ^ 
Cap) X BR(e4 ^ Wujj^)). The yellow contours are the upper bounds on (in pb) implied 
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by the H WW limits and are controlled by the dependence of our signal acceptances 
(for the WW and H WW analyses) on the H and 64 masses. The red dashed contours 
are labelled with the value of BK{H Wfiv^) cot^ (3 that leads to = 1 pb. 

Focusing on the e/i case (for which there is a larger statistics), we see that in the bulk 
of the parameter space we consider the maximum allowed WW effective cross sections are 
smaller than 10 pb and well within the allowed 2a experimental ranges (see eq. (5.4)). 
This is in contrast to what happens for H z/4Z/^ as one can see from figure 1 of ref. [1] 
where H WW constraints allow for very large effective WW cross sections in most of 
the parameter space. 

This feature is due to the different behavior of the ratio of acceptances Mnp 
the H z/4Z/^ and H 64/4 channels. This ratio controls the upper limit on the effective 
WW cross section (as we explain in appendix A of ref. [1], the larger the ratio, the larger 
the allowed cross section). Both channels have similar /^np at moderately large 
rriH and small this implies that the 10 pb yellow contours for the P4 and 64 cases are 

close to each other. As we explain below, when moving to smaller mn and larger 
the P4 ratio increases while the 64 one decreases. Because of this behavior, in the bulk of 
the parameter space in which we are interested (smaller ttih and larger we find 

large allowed a^^ values for the H z/4P^ channel but not for the H one. 

The behavior of the acceptances ratio is essentially controlled by the difference ttih — 
This difference determines the transverse mass ttit in the P4 case and the dilepton 
invariant mass mii in the 64 one. The H WW acceptance decreases for channels with 
lower rriT and increases for channels with lower mu (because the CMS Higgs cuts include 
a range for and an upper bound on mu)> The WW acceptance, on the other hand, 
is controlled by a mu > 10(15) GeV cut (for e/i and /i/i final states): for the 64 case it 
decreases at low m/f — me4, while, for the P4 one, the dilepton invariant mass is controlled by 
771^^4 and tends to always pass the cut implying a mild dependence of the WW acceptance 
on the choice of masses. In conclusion, small mn — '^^'1/4,64 implies a small acceptances ratio 
for 64 and a large one for P4. 

The discussion of the /i/i mode is similar. The main differences are that the experi¬ 
mental H WW cuts are much tighter in order to suppress Drell-Yan backgrounds and 
that the effective cross section is enhanced by a combinatorial factor of 2 with respect to 
the e/i case (see eq. (5.1)). Note that in order to obtain similar effective cross sections for 
the e/i and /i/i modes one needs to include a second vectorlike lepton family as discussed 
in section 4 of ref. [1]. 

Since the WW effective cross sections that we find in the H e4/i channel are typically 
smaller than 10 pb, we refrain from performing a detailed scan that includes mixing in the 
charged lepton sector. 

7 Contributions from SM-like Higgs boson 

In this section we consider exotic decays of the SM-like Higgs into vectorlike leptons. 
We begin by considering the h 7/47/^ ^ Will'll process. In figure 10 we show the 
effective cross sections a^^ as a function of 777^^4 G [95,125] GeV for the e/i and jifi 
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Figure 9. The contours of for the process H e^/j. W with BR(i7 ^ W{iv^) = 

BR(i7 ^ e 4 /i) X BR(e 4 ^ Wv^) = 1 (this definition applies only here) for the e/i final state (left) 
and /i/i final state (right). 
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Figure 10. Contributions of the 125 GeV Higgs decay h to the effective cross section 

as a function of rriy^ G [95,125] GeV for the e/i and modes. We set BR(zi 4 ^ W/i) = 1 and 
consider two representative values of the flavor violating Yukawa couplings = 0.02 and 0.03. 

The thick solid red line in figure 10 is the 95% C.L. upper bound from the SM h WW ATLAS 
search [39]. 


modes. Here we set BR(z/4 — 1 and consider two representative values of the flavor 

violating Yukawa couplings = 0.02 and 0.03. These values are close to the largest 

possible as one can see from the parameter scan presented in figure 11 where we show the 
— BR(z /4 W/i) plane for rriy^ — 120 GeV. The thick solid red line in figure 10 is 

the 95% C.L. upper bound from the SM h WW ATLAS search [39]. 

We see that, for the e/i mode, Higgs searches are not constraining while in the /i/i 
mode they require > 105 GeV. In both cases, the effective WW cross section cannot 
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Figure 11. Result of a parameter scan in the - BR(z /4 Wfi) plane for rriy^ = 120 GeV. 

Similar results are obtained for other choices of < 125 GeV. See the caption in figure 3 for 

further details on the scan. 


exceed 2 — 3 pb. These cross sections are far from the ranges allowed by ATLAS (blue 


shaded region) and the CMS upper bound (purple line) for the e/i final state but are close 
to the allowed ATLAS region for the fifi case (see eq. (5.4) and the related discussion). 

We do not discuss in detail the h W process because we found that it 

does not lead to appreciably large effective cross sections (typically smaller than Ipb) and 
it is severely constrained by /i ^ WW searches. 

8 Contributions from Drell-Yan production of vectorlike leptons 

In this section we discuss contributions to the effective cross section that stem from 

the following vectorlike lepton Drell-Yan production processes {I = e, ji): 


pp ^ (7, Z) ^ efef ^ ^ 2Uiy , 

pp ^ Z ^ ^ Wfih'^ . 


(8.1) 

( 8 . 2 ) 


Note that there are many more processes (involving up to four light leptons in the final 
state) that one can consider and the two modes we consider in eqs. (8.1) and (8.2) are the 
two most promising ones.^ 

The 6464 pair production channel is flavor diagonal and the Z — 64 — 64 coupling can 
be as large as the corresponding Z — I — I SM one. On the other hand, in the channel 


(8.2) the production of a single z/4 is constrained by the values of the off-diagonal coupling 
allowed by EW precision data. To quantify this effect we define 



(8.3) 


- ^2/(4cos2 0vv) 


^If more than two light charged leptons are present, the third hardest lepton must have pr < 7 GeV in 
order to avoid detection and this requirement suppresses the acceptance. 
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Figure 12. Allowed values of and for rriy^ = 110 GeV. We see that < 1.5 x 10 ^ 

and I <0.02. Similar bounds are found for different z/4 masses. 




Figure 13. The effective cross section [pb] for Drell-Yan processes. In the left panel we 

consider the channel pp ^ 4^4 2 iAiy assuming SM-like strength of 

the Z — 64 — 64 vertex, BR(e4 ^ ITz^^) = 1 and rue^^ = 105 — 250 GeV. In the right panel we show 
pp ^ Z ^ U 4 U^ £p 2 u for Ry^y^ • BR(z/4 Wp) = 10“^ and niy^ = 95 — 250 GeV. 


which shows how the production of z/4Z/^ through Z boson is suppressed compared to the 
SM process pp ^ Z ^ In Fig. 12 we see that Ry^y^ can be at most 1.5 x 10“^ for 

rriy^ = 110 GeV. 

In figure 13 we show the most optimistic values of effective cross sections for the 
processes in eqs. (8.1) and (8.2). For the latter case we set Ry^y^ • BR(z/4 -^Wp) — 10“^. 
Consequently one can see that allowed values of for the pp (7, Z) 6464 and 

pp ^ Z ^ 1 ^ 41 ^fj, channels are at most of order 1 pb and 0.1 pb respectively. 
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9 Conclusions 


We studied decay modes of a heavy CP even Higgs boson, H and H e^/j. followed 

by z /4 ^ VP/i and 64 ^ where 64 and z /4 are the lightest charged and neutral mass 

eigenstates originating from vectorlike pairs of SU(2) doublet and singlet new leptons. We 
showed that, with Yukawa couplings as in two Higgs doublet model type-H, these decay 
modes, when kinematically open, can be large or even dominant. After imposing all the 
experimental constraints, the H ^ decay channel can have branching ratio of up to 
about 35%. 

As we discussed in sections 4 and 5, electroweak precision data impose very strong 
bounds on various gauge and Yukawa couplings: the new flavor violating gauge couplings 
and have to be smaller than 0(10“^), the couplings of SM gauge bosons to 

the second family of leptons, ^.nd can deviate from their SM values by less 

than ^ 0.1%, and the flavor violating Yukawa couplings and are constrained to 
be smaller than ^ 0.05 and ^ 0.17, respectively. 

Focusing on pp ^ H ^ z/ 4 Z/^ ^ WiiU/j, we studied possible effects of this process 
on the measurements of pp WW and H WW. Contributions from this process to 
2 i2iy final states can be very large since only one W has to decay to leptons unlike in 
the case of pp WW and H WW. We present predictions of the model in terms 

of effective cross sections for pp WW and H WW in pe2u 2p2u final states from 
the region of the parameter space that satisfies all available constraints including precision 
electroweak observables and constraints from pair production of vectorlike leptons. Parts 
of the parameter space are already excluded by these measurements and thus possible 
contributions to these processes can be as large as current experimental limits. Large 
contributions, close to current limits, favor small tan [3 region of the parameter space. 

In addition, we studied correlation of the contributions to pp WW and H WW. 
We showed that, as a result of adopted cuts in experimental analyses, the contribution 
to pp WW can be more than an order of magnitude larger than the contribution to 
H WW. Thus more precise measurement of pp WW in future will significantly 
constrain the parameter space of the model. 

Furthermore, we also considered possible contributions to pp WW from H 
VF/iz/^, from similar processes involving SM-like Higgs boson and from pair production of 
vectorlike leptons. These however lead to much smaller contribution to the effective cross 
section for pp WW while satisfying limits from H WW and h WW in first two 
cases. In the case of pair production of vectorlike leptons, the cross sections are very small, 
and the contribution to the effective pp WW is at most of order Ipb. 

Finally, as we discussed at the end of section 5, the next LHC run at 13 TeV with 
100 fb“^ of integrated luminosity will be able to explore most of the parameter space 
currently allowed by electroweak precision data and H WW constraints. 
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